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We study the radiative cooling of polar and metallic nanoparticles immersed in a thermal bath 
close to a partially reflecting surface. The dynamics of relaxation is investigated at different distances 
from the surface, i.e., in the near-field and far-field zones. We demonstrate the existence of an 
oscillating behavior for the thermal relaxation time with respect to the separation distance from the 
surface, an analog of Friedel oscillations in Fermi liquids. 
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I. INTRODUCTION 

It is well known since the pioneering works of Drex- 
hage et al. [l| and Chance et al. on the molecular flu- 
orescence emission that the radiative lifetime of excited 
atoms or molecules is not an intrinsic property of matter, 
but it depends on their close environment 0, Q • This de- 
excitation process is accompanied by spontaneous light 
emission, that can be investigated by either a classical or 
quantum approach. From a classical point of view, the 
atom/molecule is considered as a simple dipole that radi- 
ates as an antenna. The radiated field is then found with 
the help of Maxwell's equations and the spontaneous de- 
cay rate is proportional to the partial local density of 
states 0,0]. On the other hand, from a quantum point 
of view, the process is described as a transition between 
different discrete states of the atom/molecule. The radia- 
tive lifetime or its emission rate is then straightforwardly 
deduced by application of Fermi's golden rule 0, 13, 01 ■ 

The radiative cooling of a nano-object is, in a way, a 
continuous version of the latter. When such an object 
is thermally excited at a given temperature, a contin- 
uum of modes centered around the hot body's thermal 
frequency is excited according to the Bose-Einstein dis- 
tribution function. As a result, these modes de-excite 
throughout heat exchanges with the surrounding envi- 
ronment. As the local temperature decreases during the 
thermal relaxation process new modes are excited (at 
smaller frequencies), changing the cooling spectrum as 
the system is driven towards thermal equilibrium. Un- 
derstanding the underlying mechanisms for such a non- 
steady dissipative process is of major importance for con- 
troling the heat flow in nanoscale systems and might be of 
practical relevance in potential applications such as near- 
field thermo-photovoltaics 3 and thermal management 
in microelectronics 

In this work, we describe the dynamic process of ther- 
mal relaxation for polar and metallic nanoparticles facing 
a substrate at a given distance. For the polar material we 
choose SiC because it supports surface phonon polaritons 
in the near-infrared. This is due to a negative real part 



of the permittivity for frequencies uj in the reststrahlen 
region, i.e. for wto < < wlO; where wto and wlo 
are the frequencies of the transversal and longitudinal 
optical phonon, respectively. As we will discuss below 
the main channel for thermal relaxation will be due to 
these surface modes. For the metal case we choose Au 
which can be described by a simple Drude model in the 
infrared region. In particular, metals like Au have a neg- 
ative real part of the permittivity for all frequencies uj 
smaller than the plasma frequency Wp which is typically 
one or two orders of magnitudes larger than wto- In pa-r- 
ticular, for Au the surface plasmon polariton resonance 
lies in the ultraviolet region so that in this case the sur- 
face modes play a minor role for thermal relaxation. In 
fact, for metals the main channel for thermal relaxation 
is due to the induction of eddy currents. 

We compare the dynamic process for both materials 
considering only particles smaller than the thermal wave- 
length Ath = he/ {k^T), where c is the vacuum velocity 
of light, 2nh is Planck's constant, fee is Boltzmann's con- 
stant, and T is the particle's temperature. Then we can 
use the dipole model based on Rytov's fluctuational elec- 
trodynamics [l^ in order to describe the heat loss of the 
particle due to radiation. This model has for example 
been used to determine the heat flux between a spherical 
nanoparticle and a flat surface (Tll - [l6| , between a spher- 
ical particle and a structured or rough surface |17H2g |. 
between ellipsoidal particles and a flat or structured sur- 
faces [21I [22I as well as between two or more spherical 
particles j23t - l26| or particles of arbitrary shape |27l . [28| . 
By using the dipole-model, we can introduce the thermal 
lifetime of a heated particle above a surface analogous 
to the lifetime of an excited atom or molecule above a 
surface. We investigate in particular the asymptotic be- 
havior of cooling rates both in near- and far-fleld regimes 
and show how the cooling depends on material properties 
of the particle and the substrate. 

The paper is organized as follows: In section II we 
introduce the dipole model for the heat flux between a 
particle and a flat surface. In section III we deflne the 
thermal relaxation time for small temperature gradients 
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and discuss it for different distance regimes for polar and 
metallic nanoparticles. Finally, in section IV we define a 
general thermal relaxation time and solve the dynamical 
problem of thermal relaxation of a nanoparticle above a 
flat surface, numerically. 



II. HEAT FLUX BETWEEN A PARTICLE AND 
A FLAT SURFACE 

Let us consider the situation depicted in Fig. [1] A 
spherical nanoparticle with a radius R is placed in front 
of a substrate at a fixed distance d. Assuming isotropic, 
homogeneous and non-magnetic materials, the nanopar- 
ticle and the substrate can be characterized by the per- 
mittivities ep and es. Furthermore, we assume that the 
particle's temperature Tp is higher than the tempera- 
ture of its surroundings Ts. Then the heat flux from 
the nanoparticle to the substrate due to the imposed 
temperature gradient can be described within a dipole 
model [TT| - [l^ if the particle radius is much smaller than 
the thermal wavelength Ath,p- Including the magnetic 
response due to the induction of eddy currents in the 
object, which is important for metals in the infrared re- 
gion m, [igI [23l . [29| . the net power dissipated by the 
dipole reads [l3| 




Figure 1: Sketch of a nanoparticle with radius R above a 
flat surface at a flxed distance d. The nanoparticle has the 
temperature Tp which is assumed to be larger than the tem- 
perature of the surrounding Ts. 



Substituting the known Green's functions for a halfs- 
pace geometrv l32l . |33| into the definition of the LDOS in 
Eq. © yields [sj 
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is the electric (i=E) and magnetic (i=H) local density 
of states (LDOS) at a distance d from the surface [30, 
HH, and ([}*(aj;z,z') is the dyadic Green tensor of the 
system (see Appendix |A|). We also have 
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where we have introduced the usual Fresnel coefficients 



which is the difference of the mean energy of harmonic os- 
cillators at temperatures Tp and Ts, respectively. Im(a^) 
and Im(a^) are the imaginary parts of the electric and 
magnetic polarizability of the nanoparticle and will be 
specified in Eqs. ([8])) and The physical interpretaion 
of expression ([T]) is the following one: The power radiated 
by a nanoemitter radiates into its environment is a func- 
tion of the LDOS of the electromagnetic field at the emit- 
ter position, which defines all possible channels for radia- 
tive and non-radiative exchanges with the background. 
Note that the above expression does not take into account 
the interaction with the image dipole (fj] . This is reason- 
able for distances d^ > i?3|es-l||ep-l|/(4|ep-H2||es + l|) 
- a condition that is automatically fulfilled when d ^ R, 
where the dipole model is valid. 
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with the z-component of the wave vector in vacuum 70 = 
•\/fcg — , the z-component of the wave vector inside the 
substrate 7s = -^/fc^es^^T?, and the wave number in 
vacuum ko — uj/c. Note, that the electric LDOS can 
be retrieved from the magnetic LDOS by interchanging 
Ts ^ fp and vice versa. Furthermore, the electric and 
magnetic LDOS contain the contribution of propagating 
modes for lateral wave vectors k < ko for which Im(7o) = 
and the contribution of evanescent modes for which 
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K > ko and therefore Re(7o) = 0. Due to tfie rotational 
invariance with respect to rotations around the z-axis, 
we can reduce ([U and (O to one dimensional integrals 
by introducing cylindrical coordinates. 

Finally, we need the polarizabilities of the nanoparticle 
to evaluate expression ([T|). In the case of a spherical par- 
ticle of radius R the polarizabilities can be derived from 
Mie scattering theory |35i]. Denoting the particle's rela- 
tive permittivity by ep and introducing the dimensionless 
variables x = fcgi? and y = ^/epk^R, one finds [1^ Issj 
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for X ^ 1, implying that the particle's radius should 
be small compared to the dominant thermal wavelength, 
R <C Ath- If we demand further that \y\ <^ 1, i.e. that the 
radius should be smaller than the skin depth at thermal 
frequencies, the above expressions reduce to 
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We have checked that with the previous expressions we 
retrieve the dipole contribution to the heat flux for an 
isolated nanoparticle, i.e., rp = rg = 0, derived by means 
of Rytov's fluctuational electrodynamics [s^ . Note that 
in the whole model we have for convenience neglected 
any nonlocal effects both in the substrate [H, i37ri4U | and 
the particle [4ll,|4i|. 



III. THERMAL RELAXATION TIME OF A 
SMALL PARTICLE 



In the absence of phase transitions, the time evolution 
of the temperature field of a nanoparticle of heat capacity 
Cp, mass density p and of volume V = AttR^/S is given 



by 
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which is just a heat-balance equation. It is important to 
point out that the substitution of ^ in the r.h.s of (fT2|) 
is apparently inconsistent, as we are using an expression 
obtained from the fluctuation dissipation theorem - and 
therefore assuming local thermal equilibrium - to calcu- 
late the cooling rate of a given object, which is clearly 
a process out-of-cquilibrium (even locally). Further in- 
vestigation, however, shows that here this does not cause 
any trouble: the characteristic time Tph for temperature 
homogenization throughout the body (Tph < 10~®s) is, 
as we shall see, much smaller than the thermal relax- 
ation time (to be defined below) . That essentially means 



that the particle cools down in a quasi-stationary fashion, 
which in turn justifies the use of H]) into (fT2|) . 

Considering a slightly heated nanoparticle with a tem- 
perature Tp = Tg+AT, a thermal relaxation time Ti (i.e., 
the inverse of cooling rate) can be defined by lineariza- 
tion of with respect to the temperature difference 
AT. The resulting linearized equation 
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In Fig. [2] we show the evolution of the thermal relax- 
ation time (TRT) with respect to the separation distance 
d from the surface for different materials obtained us- 
ing Eq. (Hi]) together with Eqs. (g])-®. We first note 
the presence of two distinct regions for the TRT. When 
the separation distance is much smaller than the thermal 
wavelength Ath,p we see that the cooling process is faster 
than for an isolated particle and it decreases drastically 
as the separation distance is reduced. On the other hand, 
in the intermediate regime for distances d ~ Ath.p we see 
that the TRT oscillates with respect to d and approaches 
for large d the (linearized) value tqj of an isolated parti- 
cle, given by [H, l36j 
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This behavior of the cooling rate demonstrates that we 
can either reduce or increase the power exchange with 
the surface for d ~ Ath.p- In order to get more physi- 
cal insight, we calculate below the asymptotic behavior 
of the power exchanged between the substrate and the 
nanoparticle in these two regimes. 



A. Cooling in the near-field regime 

The asymptotic expressions for the LDOS in the qua- 
sistatic regime, i.e., for distances z <^ Ath.p can be de- 
rived by considering only contributions for wave vectors 
K ^ fco. Then one gets the known expressions [30| 
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Figure 2: Semi-logarithmic plot of the distance dependence of 
the TRT of a nanoparticle above a substrate with temperature 
Ts = 300 K (a) for a gold nanoparticle above a gold surface, 
(b) a SiC nanoparticle above a SiC surface. The horizontal 
lines mark the asymptotic values given by Eq. (|15p . For ep /es 
we use for Au the Drude model and for SiC a Drude-Lorentz 
model m, nil . For p we use p*" and p^''^ and for the heat 
capacity we use a simple Debye model (see appendix [Cj so 
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From the structure of the previous expression we see 
clearly that the main contributions to 1/ti come from 
the particle's electric and magnetic resonances [which are 
located at ep —2 and ep — ^ oo as can be seen in Eqs. 
(|TO| - (jlll) ] and from the evanescent modes in the surface 
(at es ~ — 1). In other words, the decrease in the TRT 
can be attributed to the enhanced radiative heat trans- 
fer due to the coupling of the particle's (induced) dipole 
moment to the evanescent modes of the surface [ij, |4J] . 
In addition, it is also seen from (|18p that, depending on 
the strengths of the electric and magnetic polarizabilities, 
the lifetime can be proportional to d^ or d. For metals we 



typically observe the latter, as the induction of Foucault's 
currents gives rise to a big magnetic contribution to the 
heat transfer [2l|, [2^ , whereas for polar materials like 
SiC the electric LDOS dominates the heat flux [l^, , 
producing a. 1/d^ behavior. 



B. Cooling in the intermediate and far-field 

regimes 

On the other hand, for distances d > Xth.p wc can de- 
rive the appropriate expressions by using the stationary 
phase method (see appendix [B)) . Wc then obtain 
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We see that the presence of the substrate makes the elec- 
tric and magnetic LDOS oscillate around one half of the 
free space value, given by us^ jl-n^c? . This behavior is 
analogous to the fermionic density oscillation close to a 
defect in a Fermi fluid or to the electric charge screening 
in a pool of ions [s^l- From (fT9l) - (|20)) we get a simplified 
expression expression for the TRT, namely 
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In the case where wc have narrow dipolar resonances we 
can write 
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where /(w) is a smooth and arbitrary function that varies 
slowly on the scale of the widths of the resonances of 
Im(a^'^). Here wf are the electric (magnetic) 

dipole resonances of the particle and TVe (^h) is their 
total number, and therefore we can simplify t^^ even 
further, 



F(c.E,rs) 



+ i^(c^",Ts 



1 
2 

Im 



Im 



2ifcEd n 



1 - \/es c 

1 + ^ 4fcEcf 



dw Im(Q!^ 



1 - e 

1 -I- 4fcHd 



da;Im(a"),(23) 



where 



J^(w,rs) = 



1 / 2^3 \ d0(w,T) 



pc^y 



dT 



(24) 



Ts 



5 



ks = oj^/c, ku = i^^/c, and we assumed that only one 
resonance of each type (E and H) is present. Expres- 
sion then shows clearly the behavior of the TRT for 
large distances (and narrow resonances): it oscillates in 
a superposition of two periods while decaying as 1 /d. In 
particular, when one polarizability strongly dominates 
the other [say, F{uj^,Ts) > F{uj^,Ts)], we should see 
a clean "monochromatic" sinc-like oscillation of period 
2-K/kE- Since this is precisely the case for SiC, that is ex- 
actly what we observe in Figure 2b. That we don't have 
a similar behavior for gold is explained by (i) the fact 
that for Au particles we have iu^ > Wth for Tg « 
meaning that the electric dipole gives no contribution to 
the heat flux (because F(a;^,300) « 0), and (ii) the very 
large width of the magnetic resonance, which invalidates 
approximation (j22p and in practice smears out the mag- 
netic contribution to the TRT, resulting in the behavior 
seen in Figure 2a. 
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C. Temperature dependence 



Finally, we study the temperature dependence of the 
TRT. In Fig. [3] we show some numerical results of the 
TRT over the temperature for d = 10 fim and d = 
100 nm, i.e., in the far- field and near- field regimes ne- 
glecting the temperature dependence of the material pa- 
rameters [i^. First of all, one can observe that the 
TRT decreases with increasing temperatures Tg. Hence, 
the hotter the environment, the faster the particle cools 
down (reminding that, in the linear approximation, Tp = 
Ts + AT with AT < Ts) . Note that in the given temper- 
ature range of 50 — 500 K the TRT varies enormously. For 
Au (SiC) particles at d = 10 /im it can be as big as 4 hours 
(98 minutes) at Tg = 50K, dropping all the way down to 
around 5 seconds (266 milliseconds) at Ts ~ 500A". Even 
more interesting is what happens in the near-field regime 
(d = 100 nm), where the TRT variation is not so impres- 
sive for gold (from about 240 milliseconds at Ts = 50A' 
to 58 milliseconds at Ts = 500K) but still very large for 
SiC (around 5 minutes at Ts = 50K and 0.3 millisec- 
onds at Ts = 500K). This difference at close distances is 
readily understood in terms of the coupling to the surface 
modes of the substrate. For SiC we have surface waves 
that lie in the frequency region close to the maximum of 
the thermal function F{uj, 300) defined in ([M]) . and so the 
heat exchange is quite large at those temperatures. By 
decreasing Ts, however, we shift T(w,Ts) toward lower 
frequencies and strongly suppress the participation of the 
surface modes in the heat transfer, leading to the huge 
increase in the TRT. By contrast, the surface modes for 
Au lie at much higher frequencies, which means that they 
don't play any significant role in the heat fiux in the con- 
sidered range and therefore can't affect much the TRT. 
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Figure 3: Semi-logarithmic plot of the temperature depen- 
dence of the TRT of a nanoparticle above a substrate at a fixed 
distance d = 10 fim and d — 100 nm for (a) a gold nanoparti- 
cle above a gold surface, (b) a SiC nanoparticle above a SiC 
surface. The particle radius is in both cases R — 10 nm. 



IV. DYNAMICS OF RELAXATION 



The TRT defined in Eq. ^ is only useful for small 
temperature differences AT, since it is based on a lin- 
earization procedure. Here, we consider now the full dy- 
namical process regarding the cooling of a particle by 
solving numerically the nonlinear Eq. (jl2p . Here, we 
need to define the TRT in the general case: Let us as- 
sume that the nanoparticle is heated up with respect to 
its surrounding (Tp > Ts), then we define the thermal 
relaxation time of the particle as the elapsed time r the 
particle needs to cool down to Ts + AT/e. This defini- 
tion ensures that for very small temperature differences 
we retrieve the linearized relaxation time ti defined in 

Eq. mi). 



A. Cooling dynamics for an isolated particle 

For an isolated particle we already saw that = 
= Dq/2, so that the cooling dynamics is governed by 
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the following equation 
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(25) 

In Fig. m we show the results obtained by a numerical in- 
tegration of Eq. (|25|l using a Runge-Kutta method with 
adaptive time steps. As for the initial conditions, we 
consider two different cases: (i) The solid lines show the 
TRT (as defined above) for Ts = SOOi^ and Tj.(t = 0) = 
SOOif + AT, while (ii) the dotted lines correspond to 
Tg = 300A' - AT and Tp(t = 0) = SOOif. It can be seen 
in Fig. m that in the first case the TRT decreases with 
respect to AT, whereas in the second case it increases. 
This somewhat unexpected result shows that the cooling 
rate strongly depends not only on the temperature differ- 
ence between particle and environment, but also on the 
absolute temperature of the environment. The first sce- 
nario can be qualitatively understood from the fact that 
by heating up the particle (while Ts is kept constant) 
the radiated power increases nonlinearly with Tp, lead- 
ing therefore to a faster cooling. In the second case the 
radiative power of the particle at < = is kept constant 
while Tg is cooled down. Thus the increasing difference 
between Tp and Tg leads to an increase of the time the 
particle needs to cool down to Tg -I- AT/e. 



B. Cooling dynamics of a dipole near a surface 

Now we look at the evolution of the temperature of 
a slightly heated nanoparticle near a surface. In Fig. [5] 
we show the time evolution of Tp for a gold nanoparticle 
with R = 100 nm for a distance d — 500 nm above a gold 
surface and a SiC nanoparticle above a SiC surface. As 
previously we consider a fixed particle temperature or a 
fixed surrounding temperature. As for the isolated parti- 
cle, the exact relaxation times are in very good agreement 
with their linearized counterparts for AT — )- 0. We ob- 
serve that due to the presence of the interface the TRT 
is much smaller than for the isolated particles. Apart 
from this we find the same qualitative behaviour as in 
the isolated case in Fig. H) 



V. CONCLUSIONS 

We have investigated the thermal relaxation of a 
nanoparticle close to an interface. In particular, we have 
introduced a relaxation time for the regime of small 
temperature differences between the particle and the sur- 
face and a generalized relaxation time t for arbitrary 
temperature differences. Equipped with these definitions 
we have studied numerically the thermal relaxation of 
gold and SiC nanoparticles close to a gold or SiC sur- 
face, respectively. We have shown that the relaxation 
time behaves similarly to the lifetime of an atom close to 
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Figure 4: Semi-logarithmic plot of the TRT of a gold (a) and 
a SiC (b) nanoparticle with radius R — 100 nm in free space, 
as a function of the temperature difference AT = rp(0) — 
Ts (see text for further details). The horizontal line is the 
corresponding linearized TRT Tg. 



a surface, varying rapidly for distances smaller than the 
thermal wavelength and finally droping to zero. In addi- 
tion, we found that the thermal relaxation time ti is very 
sensitive to the temperature and can vary drastically - by 
several orders of magnitude - in the temperature range 
of 50 — 500 K. We have identified all mechanisms which 
govern the thermal relaxation from subwavelength sepa- 
ration distances (in near-field regime) all the way to long 
separation distances (in far field regime). 

This work can be generalized to the heating/cooling of 
a nano-object immersed in a more complex environment 
and exploited for the thermal management at nanoscale 
of complex plasmonic systems. Finally, a time-resolved 
scanning near-field optical microscopy as the one used in 
Refs. 1131 could be used to measure the time evolution 
of a substrate temperature in non-equilibrium situation 
by detecting, at the scale of TRT, the time variation of 
diffracted signal. 
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Appendix A: Green's dyadic for a flat interface 

The Green's dyadic function with observation point 
r — {x,y,zy and source point r' = {x',y',z'y above 
the flat surface, i.e., for < z < z' can be expressed 

as mill H 



G(°)(r,r';w) = 



d^K iei'^-(^-''') 

(2^ 2^0 



Oftn 



(Al) 



where is the unit vector in z-dircction and (g) symbol- 
izes the dyadic product. The tensors 1 and K, are defined 



where 



3 = {s.p} 
j={s,p} 



K 



1 



2\t 



^P^**) = —77(^2:70, fcy 70, ) 



(A2) 
(A3) 

(A4) 
(A5) 
(A6) 



are the polarization vectors for s- and p-polarization. 
Note that these vectors arc always orthogonal, but only 
have a unitary norm for propagating modes with k < feg. 
The reflection coefficients Vg and rp are the usual Fresnel 
coefficients in Eqs. ([6|) and d?]). 



Appendix B: Density of states in the intermediate 
and far-field regimes 

We describe here the asymptotic expansion of the TRT 
in the far-field regime. To this end we start from the 
expressions for the electric and magnetic LDOS from 
Eqs. and (O and take only the propagating modes 
into account. The distance independent part can be sep- 
arated from the distance dependent part, so that the 
LDOS can be written as 



Do 
2 

Da 
2 



Re 



Re 



'° d^K Do 



(27r) 4fco7o 



d^K Do 



fr2 



A-2 



(27r) 4fco7o 



2k' 



k' 



k' 



(Bl) 
(B2) 



where Dq = a;^/(7r^c^). That means for the distance 
dependent part we have two-dimensional integrals of the 
form (i = E, H) 



h = Re / d''KFi{K)e''''^'^> 



(B3) 



which have a highly oscillating integrand for 4'{k) larger 
than one, which corresponds in our case to distances 
larger than the wavelength. By using the stationary 
phase method [iot such integrals can be approximated 
by 



Re- 



2n 



i\/ ac — b'^ 
where Ks is defined such that 



d(j) 
dk„ 



= 



(B4) 



(B5) 
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and 



model |5C 



9^ 
dkl 



b = 



and c — 



(B6) 



In our case = 27od so that Kg ~ (0, 0)* and a = c = 
2dlkQ and 6 = 0. Hence, by inserting these results into 
Eq. (jB4[) wc find the far-field expressions in Eqs. (flQll and 

(laa. 



Appendix C: Heat Capacity 

In order to describe the heat capacity of the two con- 
sidered materials (Au and SiC) we use a standard Debye 



NAkj 
' M 



T 



eo/T 



(e^ - 1)2 



(CI) 



where TVa = 6.022 • lO^'^moP^ is Avogadro's constant, 
M is the molar mass and ria the number of atoms or ions 
per unit cell, and 0d is Dcbye's temperature. For con- 
venience we assume that Cy ~ Cp. When using the pa- 



rameters 
and 



170 K, 71^" = 1, M^" = 196.97 gmol" 



HOOK, nf^ = 2, AfSiC ^ 40.1 gmol"^ we 
get values for Cp which arc in very good ag reement with 
experimental and numerical results j5ll454| . 
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